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Abstract 

We investigate the structure of Polchinski's formulation of the flow equations for 
the continuum Wilson effective action. Reinterpretations in terms of I.R. cutoff 
greens functions are given. A promising non-perturbative approximation scheme is 
derived by carefully taking the sharp cutoff limit and expanding in 'irrelevancy' of 
operators. We illustrate with two simple models of four dimensional X(f 4 theory: 
the cactus approximation, and a model incorporating the first irrelevant correction 
to the renormalized coupling. The qualitative and quantitative behaviour give 
confidence in a fuller use of this method for obtaining accurate results. 
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1. Introduction. 



It need hardly be stated that an efficient method of performing accurate ab initio 
non-perturbative continuum calculations in realistic quantum field theories would 
prove extremely useful. What is needed is to be able to produce a sequence of 
better and better approximations in the sense that they can be seen to converge. 
The higher approximations must be calculable without inhuman effort, and if the 
method is to be generally useful it should apply even if there are no obviously 
identifiable small parameters to control the approximation. This paper we hope is 
a first step towards such a scheme. 

The scheme indeed seems to work, judging by the preliminary model approximations 
presented in this paper. Since the majority of the paper discusses more formal 
aspects we limit ourselves here to these simple models, but the apologia should be 
added immediately that these are not enough to be sure that the requirements of 
convergence and calculability really are satisfied. A proper demonstration of this 
requires that the method be applied to some real non-perturbative problem, such 
as triviality bounds on the Higgs mass with and without 1 the top. This research is 
underway. One should also add that since the method is tied to momentum cutoff, 
gauge-invariant theories present special problems of their own. 

We use Wilson's exact renormalization group[l] applied to the continuum field 
theory with some physical cutoff, in the particularly simple form introduced by 
Polchinski[2]. Our reasons for using the exact renormalization group are described 
in sect. 2: it is the only framework we know of which avoids the difficult U.V. 
problems afflicting other approximation methods. In particular we contrast with 
attempts to use approximations to Dyson- Schwinger equations, and briefly compare 
other alternatives. 

In the approach of ref.[2] the results are phrased in terms of an effective action 
Sa[¥>] with an intermediate-scale momentum cutoff A; from this we can construct 
the greens functions and S matrices of ultimate physical interest but they will 
be restricted to having momenta less than A. Originating with Keller, Kopper 
and Salmhofer[3] it was noticed that the same equations as Polchinski's are also 

to compare with lattice results 
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satisfied essentially 2 by Wa[J], the generator of connected greens functions with 
I.R. momentum cutoff A. Therefore in this case as A — > one is left with greens 
functions defined for all physical momenta. 

It is immediately clear that the two constructs must be closely related. While this 
is evidently understood by workers in the field, 3 no-one seems to have precisely 
spelt out the relation. Therefore we start with this in sect. 3. The fact that, as 
we make clear, the vertices of Sa [<£>] are essentially both vertices of the Wilson 
effective action and I.R. cutoff connected greens functions, depending on whether 
the external momenta are all below or above A respectively, proves very useful both 
in our analytic investigations but also interpretationally. 

We now describe the results reported in the rest of the paper. We endeavour to give 
here, as much as possible, an intuitive explanation of the results; they are proved 
properly later. In some places the reader might find the intuitive explanations 
unilluminating, separated as they are from the details and exact definitions. If so 
the reader should ignore them for the moment and return to them later, having 
read the main body of the paper. 

Turning to the flow equations we note that they are very sensitive to the form of the 
intermediate cutoff A, even becoming ambiguous in the limit of sharp cutoff. This 
is simply because they describe the effect of integrating out momenta around the 
cutoff (p«A), where the effective vertices are changing rapidly. Therefore, as the 
equations stand, they are inappropriate for developing approximations. We resolve 
the limit by finding the general solutions, parametrized in terms of functions that 
are smoothly varying at p = A. These functions can be chosen to coincide with 
the 1PI (one particle irreducible) vertices of Ta^ ], the Legendre effective action 
with I.R. momentum cutoff. Intuitively this follows if one recalls that the vertices 
of S\ correspond to I.R. cutoff connected greens functions. These have a tree 
expansion in 1PI vertices connected by I.R. cutoff full propagators. Thus Sa can be 
reconstructed from the 1PI greens functions. But these latter are smoothly varying 
at p = A because in them all I.R. cutoff propagators are integrated over internal 
momenta, so there are no terms left to vary rapidly as the external momenta range 

2 actually the external legs must be multiplied by inverse bare propagators 

3 see for example ref.[4], ref.[5] and references therein 
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over p A. Performing the Legendre transform gives flow equations for the 1PI 
vertices. (The flow equations are also derived in refs.[8][5]. See also the "Note 
Added" at the end of the paper.) In fact we will show that the vertices of the 
Legendre effective action Ta are the same as those of the Wilson effective action S\ 
when all partial sums of external momenta are less than A (because in this range 
all the connecting I.R. cutoff propagators in the tree expansion vanish). Thus the 
dual interpretation alluded to above is reflected here too. In fact Pa[<£> c ] has an 
interpretation as a Wilsonian quantum effective action obtained by integrating out 
purely quantum modes with momenta p > A; S\ being obtained by integrating out, 
in addition, the tree level excitations of the "classical field" tp c with momenta p > A. 

One final simplification of the structure of the equations is possible before turning 
to approximations. Up until now all manipulations have been regulated by an 
overall momentum cutoff Ao and indeed this appears explicitly in the flow equations. 
One may ask however, why equations that only require momenta p ~ A should 
depend on A » A at all. In answering this question we show that there exists a 
reparametrization invariance of the flow equations that allows us to change the cutoff 
Ao without at all altering S\ [ip]. The Ao = oo case gives the maximally analytic 
resolution of the original flow equations. Solutions are obtained by setting boundary 
conditions at some finite scale A . The reparametrization is simply a tree diagram 
expansion of the 1PI vertices with overall momentum cutoff Aq in terms of the full 
propagators and 1PI vertices with a higher overall momentum cutoff. Surprising as 
this result may seem it can be understood intuitively from our previous comments: 
On the one hand the 1PI vertices with the lower overall momentum cutoff are not 
1PI with respect to the higher overall momentum cutoff, because the expansion of 
connected greens functions in terms of 1PI greens functions was performed only for 
propagators with momenta less than the lower cutoff. On the other hand we know 
that we may compute the Wilsonian effective action at cutoff A = A in terms of 
tree diagrams constructed from the Legendre effective action with the higher overall 
momentum cutoff and I.R. cutoff A . This Wilson effective action corresponds to the 
bare action for an effective theory with U.V. cutoff A . Using this to construct the 
effective Legendre effective action by integrating out quantum modes with momenta 
less than Ao will give the latter a tree expansion also. 
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Finally in sect. 4 we discuss approximations. Needless to say the equations are still 
not solvable exactly: there are an infinite set of them, they are non-linear, and they 
contain explicitly I.R. cutoffs on all propagators. It is the last property that makes 
them particularly difficult for analytic (and numerical) methods. For example even 
in the sharp cutoff limit, one-loop diagrams with external momentum dependence, 
are not doable exactly. It is therefore necessary that any approximation simplify the 
effect of these cutoffs. We consider several approximations, but the most promising 
is to make some expansion in the external momenta around p = 0. It turns out 
that a truncation coinciding with the simplest case, where all external momentum 
dependence is discarded, has been discovered and rediscovered several times[6]-[9] 
(probably only a partial list, for further references see the "Note Added" at the end 
of the paper). This truncates to an (uncontrolled) approximation for the effective 
potential. Nevertheless the results are already encouraging, incorporating as they 
do the exact one- loop answer plus some higher corrections, and give reasonable 
qualitative descriptions and numerical results for critical indices etc. in scalar field 
theory for gaussian and non-gaussian fixed points and dimensions 2 < D < 4 [7] [10]. 
The problem of sensitivity to the form of the cutoff, discussed above and in sect. 3, 
can be avoided at this level. 

Considering now the higher orders, it becomes important whether one uses a 
smooth[8][9] or sharp[6][7] cutoff. With a smooth cutoff 6 £ (q, A), momenta q in 
the flow equations are integrated out over a shell of radius A and thickness 2e. 
Evidently we must have e < A. The momentum expansion corresponds to a local 
derivative expansion in the effective lagrangian with a radius of convergence p ~ e 
(from expanding terms such as e (|q + p|, A) with q pa A). Since these expansions 
are substituted back in the flow equations where they are themselves averaged over 
p pa A, we must have e > A for convergence. The two inequalities together imply 
that the expansion method converges very weakly if at all. One also finds that 
the coefficients of the derivative terms are very sensitive functions of the form of 
the cutoff [9] , being proportional to integrals over shells q pa A of higher derivatives 
of e {q, A) (which in general must be performed numerically). This makes their 
interpretation difficult. Formally this sensitivity disappears in the limit A — > but 

4 See ref.[8][10] for an incorporation of one-loop wavefunction renormalization. 
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also in this limit the momentum expansion breaks down completely. 

For all these reasons (explained further in sect. 4) we choose instead to take the sharp 
cutoff limit. Many integrals in the momentum expansion (averages over 3-spheres 
radius q = A) are now doable exactly, and the resulting flow equation is a set of 
simultaneous first order differential equations which are easily solved numerically. 
Moreover we expect the radius of convergence to be the natural one (the minimum 
distance from the region q > A to the first singularity in complex q space) . However 
we no longer have the luxury of a local effective lagrangian. The non-localities are 
non-analytic terms in the momentum p^ and are straightforwardly related to the 
fact that both the cutoff and momenta are now precisely determined: momentum 
modes <p(q + p), where p « q = A, are integrated out if and only if p.q > 0, no 
matter how small p is. Thus the expansion is in the length p (not p 2 ) and when 
several momenta are involved the coefficients are non-trivial functions of the angles 
subtended between them. 

A sensible method of approximation is now to truncate at some order p m and at 
some n-point function (setting higher order in n and m arbitrarily to zero), the levels 
being determined by how relevant (in the colloquial sense) the higher corrections 
are to some physical quantity of interest. As a rule of thumb one might expect this 
to coincide with an expansion in increasing order of irrelevancy, as in the technical 
sense. We investigate in sect. 5 two simple models, as a first check for convergence, 
calculability and any difficulties of interpretation in this scheme. 

The first model is the analogue of the ladder (or cactus) approximation for four 
dimensional \(p 4 theory. This can be solved analytically and the flows, fixed points 
etc. understood. We note that, choosing a certain range of negative bare mass- 
squared TOq and bare coupling Ao, and working in the unstable symmetric phase, 
makes the model track towards the tachyon singularity, giving a peculiar continuum 
limit. However we expect zero radius of convergence here. The problem is cured by 
expanding about some other point <p = const. ^ e.g. the true minimum. Otherwise 
the results are all physically sensible. In particular we find a sort of "triviality of 
mass" , namely that it is impossible to have a small renormalized mass-squared m 2 
(compared to the cutoff Aq) if the coupling A is greater than a critical value A c = 
8tv 2 !=a 79. This value, interpreted as the maximum renormalized coupling, agrees 
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unreasonably well with previous precision lattice calculations [11]: A max = 78 ± 3. 

In the second model we investigate non-perturbative corrections to the coupling; 
the raison d'etre of the method. Here we make the unnecessary ansatz, purely 
to keep it technically simple, of setting all self-energies to zero. Furthermore we 
truncate by setting the 6-point function to zero. The perturbative corrections to 
the coupling incorporate the exact one-loop result and diagrams at two loops that 
are not associated with wavefunction or mass renormalization (cf. fig. 5). The 
nested two-loop diagram is not calculable exactly (because of the I.R. cutoffs) but its 
contribution to the /9-function is calculable to high accuracy as a rapidly convergent 
numerical series, by expanding the inner loop in the outer loop's momenta. This is 
a nice model confirmation our general philosophy. Ordinarily such a method would 
fail disastrously: the expansion being both increasingly I.R. and U.V. divergent in 
the inner and outer loop respectively. It works here because inner loop integrals are 
I.R. cutoff by the outer loop's momentum. Performing the momentum expansion 
to next-to-leading order on the non-perturbative equations, we construct a term 
of the form ip 2 \/n </? 2 in the effective lagrangian, giving a /3-function which is an 
asymptotic power series in A(A) incorporating the exact one-loop contribution and 
99.1% of the 2-loop term. Integrating the non-perturbative equations, they focus 
in on the gaussian fixed point: A(A) — > as A — > 0. The correction to the one- 
loop result is small. For example in the "scaling region" [11], A(A) < A s ~ 5, we 
find that the percentage correction to A(A) depends on A(A) as aA(A) where a = 
— 1.7% (—5.6%) for bare couplings Ao = 20 (1000). In fact there is a surprising 
amount of the physics of triviality captured by this model: one can show that 
the non-perturbative effects are all cutoff dependent and measure their magnitude; 
compute the RG improved Symanzik terms, the large order behaviour of the series 
and the renormalon ambiguity. Such investigations however take us too far from 
the main purpose of this paper and will be reported on elsewhere. 

Taken together these preliminary tests suggest that we have the technique under 
control, while the smallness of the corrections to the one-loop case — even for 
very large bare couplings — gives confidence in a fuller use of these approximation 
methods for obtaining accurate results. 

6 



2. Renormalizability. 

As noted in the introduction our aim is to find a method of continuum calculation 
in realistic quantum field theories consisting of a sequence of better and better 
approximations, calculable without inhuman effort, and applicable even if there are 
no obviously identifiable small parameters to control the approximation. 

Our first thought in this direction was to use the Dyson- Schwinger equations. It 
is well known that these are solvable in the large N approximation. However one 
quickly runs up against the above desiderata: while higher corrections in realistic 
cases, e.g. the Higgs sector with 0(4) — > O(iV), are calculable in principle, in practice 
the difficulties are such that the results are essentially restricted to N = 00 [12], 
while it is easy to envisage cases where the true N is definitely not large enough 
for the asymptotic expansion in 1/N to give good estimates to any order. Instead 
we tried to justify and improve the "ladder 5 approximation" and its (arguably 
better) variants [13], 6 by incorporating directly some more of the infinite set of 
Dyson-Schwinger equations. It is here that one runs into what appears to be 
an insuperable problem: In order to make progress one must truncate the higher 
equations somehow, but the results (beyond the ladder) seem always not to be 
even perturbatively renormalizable. This is because one throws away with the 
truncation certain sets of divergent diagrams, but keeps the simpler counterterm 
diagrams generated by lower terms in the Dyson-Schwinger equations. This problem 
is only seen beyond the ladder because it is only here that the coupling constant 
starts receiving divergent corrections to all orders in perturbation theory. The 
counterterms to which we refer are those generated by replacing the bare vertices in 
the lower terms of the Dyson-Schwinger equations by the required power series in 
renormalized coupling. While these cancel divergences in diagrams that are kept, 
they also generate counterterm diagrams corresponding to diagrams that have been 
thrown away with the truncation. Since the counterterms no longer have divergences 
to cancel, they contribute divergent terms (in particular some of the form In pin A 
etc. ) which destroy perturbative renormalizability. It is not clear to us that this 

a.k.a. "rainbow" or (for A</? 4 theory) "cactus" 

got by keeping only the Dyson-Schwinger equations for the self-energy and intro- 
ducing some ansatz for the vertex 



problem can at all be solved, but more importantly even if it can, the resulting 
sequence of approximations are unlikely to be easily calculable. 

(A number of authors have tried to improve the ladder ansatz by perturbatively 
correcting the kernels, e.g. the photon propagator in strong coupling QED. This only 
worsens the approximation however because the perturbative series is asymptotic 
and inapplicable to the strong coupling regime. It is easy to explicitly verify this 
with a simple model (0 + dimensional \(p 4 theory) but this would take us too far 
from the main thrust of this paper.) 

This experience teaches us an important lesson however, which is applicable to many 
a scheme involving truncations [14]: namely that perturbative renormalizability is 
not guaranteed. If it is truly necessary to test all attempts to repair or improve these 
schemes by some analogue of the detailed graphical and convergence arguments used 
in the classic proofs of perturbative renormalizability, then this is rather depressing, 
not the least because they likely fail the test and in any case non-perturbative 
renormalizability — even for very small coupling — would remain an open question. 

It is for the above reasons that we turned to the exact renormalization group as 
the appropriate framework for non-perturbative approximations: As pointed out by 
Polchinski in a beautiful paper [2], in the Wilson renormalization group framework 
the perturbative renormalizability of e.g. four dimensional A</? 4 theory, is essentially 
obvious. It is only necessary to show that the operators we know are irrelevant (i.e. 
have negative scaling dimension) at zero coupling, remain so at very small coupling. 
But this is so because the right hand side of the Wilson flow equations dS\/dA = • • •, 
where A is an intermediate cutoff for the effective action Sa, is a smooth function 
of the coupling. Indeed the detailed proof[2] only requires very simple bounds on 
the right hand side which do not involve any cancellations between terms, and thus 
it is clear, both intuitively and in detail, that truncations of the flow equations are 
also perturbatively renormalizable. Of course what happens for large couplings is a 
different matter since here renormalizability even in the full theory is brought into 
question, but at least one can be sure that this framework is sensible in the small 
coupling regime. 

Actually the problem of renormalizability in truncations has its analogue in diffi- 
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culties with any approximation, in the following sense: If we are only modelling 
the high energy behaviour approximately then the divergences of quantum field 
theory will tend to ensure that the approximation is infinitely bad in the limit that 
the cutoff is removed. It is such an unwelcome sensitivity to high energies (together 
with the inability to compute any but the simplest approximation) which is the main 
reason for the failure of variational methods in quantum field theory [15]. Lattice 
methods suffer too in the sense that nearly all the numerical effort is expended on 
computing effects of order the lattice spacing, the time required and the number 
of lattice points required diverging as the parameters are tuned to the continuum 
limit. 

The methods described in this paper appear not to suffer from this general 
problem; on the one hand the calculations are protected from U.V. divergences 
by the intermediate cutoff A and on the other, even approximations to the exact 
renormalization group focus the effective lagrangian — erasing errors in irrelevant 
degrees of freedom and automatically absorbing them in the low energy definition 
of relevant and marginal couplings. 

3. The Flow Equations and Their Structure. 

Throughout the paper we work with a one-component scalar field <£>(x): the 
generalisation to more components and fermions is straightforward. As usual in 
momentum space, we factor out and evaluate the momentum conserving 5-functions 
so n-point greens functions are written as ~ G(p±, • • • , p n ) and are defined only 
when pi + • • • + Pn = 0, while in addition in two-point functions we solve pi = 
— P2 = p and recognise that they are functions only of p = |p| e.g. the propagator is 
written A(p). We use the condensed notation wherever convenient: Lp.J = ip x J x = 
Jd D x v?(x)J(x) where D is the dimension of euclidean space-time. Similarly the 
propagator A(x, y), and other greens functions of two arguments, will be regarded 
as a matrix so cp.A~\(p means (p x (A~ 1 ) xy (p y = Jd D p/(2n) D (p(p) [A(p)] _1 ip{— p). 
The partition function is assumed to be regulated by an overall momentum cutoff 
Ao; this will be made explicit in a moment. 

Let us start by observing that the partition function for the scalar field ip with 
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propagator A and arbitrary bare interaction S\ Q [<p] , 

Z[J} = Jvtp exp{-^<p.A-\(p-S Ao [(p] + J.<p} (3.1) 

can be rewritten in terms of two propagators and two fields as 

Z[J} = JV(p > Vcp < exp{-|^>.A;V> - 5^<-A<V< 

- Sa [v?> + <P<] + J-(<P> + <P<)} 
(up to a multiplicative factor which we always ignore), where 

A = A< + A> and <p = ip > + <p < . (3-3) 

This is true whatever 'partition' of A is used and is obvious in perturbation theory, 
since every propagator in each graph is just repeated twice — once with A< and 
once with A>, but it is also true non-perturbatively. (Substitute <p > = <p — <p < 
followed by (p < = (p' K + (A</A).<p, and integrate out the now gaussian <//<). This 
trick was used in ref.[16] for a nice constructive proof of Polchinski's flow equations. 
Here we will use it to further investigate the properties of the effective action. 

Write 

A>(p) = [0 e (p, A) - 9 £ (p, A )] A(p) and A<(p) = [1 - 6 e (p, A)] A(p) , 

(3.4) 

where 9 e (p, A) is a smooth cutoff function, bounded above (below) by one (zero), 
satisfying 9 £ (p, A) for p < A — e and ^ e (p, A) ps 1 for p > A + e. Thus A> (p) is the 
propagator cut off from below by A and above by Ao, while A<(p) is the propagator 
cut off from above by A. We have taken the opportunity here to make explicit the 
overall momentum cutoff by replacing A(p) in (3.1)-(3.3) by [1 — 9 £ (p, Ao)] A(p). 
The intention eventually is take the sharp cutoff limit: 

9 e (p, A) ^9{p- A) as e^O (3.5) 

and thus to identify, by (3.2), </ 7 <(p) with the low momentum modes p < A and 
</?>(p) with high momentum modes p > A. But we keep A<(p) and A>(p) non-zero 
for all p at intermediate stages, to inject an element of rigour into the approach. 
(Observe that for limit (3.5) the inverses in (3.2) are not defined). 
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Integrating out the nascent high momentum modes <p > in (3.2) we have 

Z[J]=jVi P< exp{-±<p < .A- 1 .<p < -S A [A > .J + <p < ] + J.<p < + lJ.A > .J} (3.6) 

for some functional S A . To see this, isolate in (3.2) the integral over (p > dependent 
factors, and substitute <p > = <p — </?<: 

Z[J] = J^Pk exp{ — |(/7 < .A~ 1 .(/? < } Za [</?<, J] where (3.7) 

Za[v<, J] = J £V> exp{-i(^>.A; V> - Sa [<P> + ¥<} + J-{<P> + </?<)} (3.8) 
= exp{-i^<.A;V<} /^V exp{-|^.A;V-^AoM+^(^ + A;V<)}(-3.9) 

Integrating over <p gives 

Z A [ip Kl J] =exp{ij.A>.J+ J.</?<} {exp-|(J + A;V<)-A>.(J + A;V<) x 

exp-^Aot^y] expi(J + A;V<)-A>.(J + A;V<)} • 

Performing all the derivatives in S Ao [jj] and noting that the ^j's are replaced by 
either A>. J + <p < or the differential of this (A>), we have for some 5a 

Z A [v?<, J]=exp{±J.A > .J + J.<p < -S A [A > .J + <p < }} , (3.10) 

proving the assertion. 

Now consider the limit (3.5). If we also insist that J(p) = for p > A, so that J 
only couples to the low momentum modes (as in ref.[2]), we see that in (3.6) all the 
J's drop out except for </.</?<, and thus S A coincides with the Wilsonian effective 
action 7 : it is the same form as (3.1) with Ao i— > A. We can keep this interpretation 
for general J(p) if we interpret J in (3.1) as a space-time dependent one-point 
coupling, analogous to a non-constant external magnetic field. 

On the other hand if we set <f < = in (3.8), we have a standard partition function 
for a field (p > with an I.R. cutoff A imposed. Thus using (3.10) and defining 

W A [cp < ,J}=lnZ A [cp < ,J] =±J.A > .J + J.(p < -S A {A > .J + ( P< } , (3.11) 
7 in the same form as ref.[2] 
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we have that Wa[0, J] = \ J.A>. J — 5a [A>. J] is the generator of connected greens 
functions with I.R. cutoff A. This is the relation advertised in the introduction: 
we see that in the limit (3.5) the support of Sa neatly separates into low momenta 
(provided by <^<) where it is the Wilsonian effective action, and high momenta 
(provided by A>. J) where it is related to the generator of connected greens functions 
as above. To make the relation completely explicit, define $ = A>. J + <p< and 

s(pi '"'' p "' a) = «*( Pi )...«*(p„) ' (312) 

Writing the I.R. cutoff connected greens functions as 

S n W A [0, J] 



G(pi, • • • , p n ; A) 



5J(pi)---5J(p„) 
we obtain in the limit (3.5) 

A > all pi : S(pi, • • • , p n ; A) = vertex of effective action at momentum 

scale A 

n 

A < all Pi : 5(px, • • • , p n ; A) = -G(pi, • • • , p n ; A) J| A" 1 ^) for n > 2 

i=i 

S(p; A) = [A(p) - G(p; A)] at n = 2 

(3.13) 

(Of course in the second relation we need also all pi < A ). 

At this stage we interrupt the exposition with some pedantry. Firstly observe 
that, in the limit (3.5), v ? <(p) = for p > A is not imposed in (3.10), but only 
fluctuations satisfying this give non-zero contributions to (3.7). Our later proofs can 
be shortened by using this observation to note that Sa in (3.10) may therefore be 
regarded as over-parametrized. Indeed all that is needed to identify the Wilsonian 
effective action as above, and effectively to make all our later observations also, is the 
limit A>(p) ^0 (as e — > 0) for p < A. Thus (in the limit) we require A<(p) = A(p) 
for the low momentum modes p < A, while for p > A the propagator A<(p) can 
for example have a high energy 'tail', giving cp < support on all momenta. We will 
not continue to make or use such comments, because they neither clarify nor are 
they important for our present purposes. The same applies to modifications of 
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our observations to fit the appropriate inequalities and approximate relations for a 
smooth cutoff (3.4), and where the approximate relations become exact for 9 e {p, A) 
an exact partition of unity. We trust the reader interested in these cases will find 
it straightforward to modify the proofs as appropriate. 

Let us define 



K A (p) = --^-A > (p) = d £ (p,A) A(p) (3.14) 



so, by (3.5), 



S e (p, A) = -j^ e {p, A) -> S(p - A) as e^O . (3.15) 
From the definition (3.8) we have 

which, on substituting (3.10), gives 

d 1 jSS A 5S A ( S 2 S A \\ 

M SAm = 2{l¥- KA -l¥- tr { KA -s^)S ■ (3 - 17) 

(N.B. again, in this equation K A and the $ two-point function are regarded as 
matrices). This the Polchinski equation [2]. It can also be written in linear form as 

m exp - %m = 4 tr ( K * jm exp ~ s * m ) ■ 

but this doesn't seem to be useful. Expanding (3.17) using (3.12), we obtain 

As( pi) ... )Pn; A) = S(-Y> l ,I l ;A)K A (P 1 )S(P 1 ,I 2 ;A) 

{/l ' /2} (3.18) 
1 r d D q 



- 2 / T^Xd K A (q)S(q, -q, Pi, • • • , p n ; A) 



(27T) 

In here I\ and I 2 are disjoint subsets (ii fl 7 2 = 0) of the momenta such that ii U I 2 = 
{Pi; • • • > Pn}- The sum over {ii, J 2 } means a sum over all such disjoint subsets, but 
pairs are counted only once i.e. {Ii, I 2 } = {I2, h}, utilising the Bose symmetry of the 
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Fig.l. The Polchinski equation for the vertices. These are drawn as open circles, 
while the two-point function K\ is drawn as a black dot. 

equation. The momentum Pi is given by momentum conservation as Pi = J2 Pi eii P« 
. The equation is best appreciated graphically — c.f. fig.l. 

Eqn.(3.18) gives the change in the vertices of the effective action, under an 
infinitessimal lowering of the intermediate cutoff A, as an infinite set of simultaneous 
first-order differential equations in A. Evidently they can in principle be used to 
define and solve the theory, the initial boundary conditions being given by S\ = S\ 
at A = A . With just two terms, totally symmetrized, the first representing tree- 
level contractions between vertices with the same or smaller numbers of legs, and 
the second a one-loop integral over a vertex with two more legs, it is about as simple 
as one could hope for. 

They are not however simple from the point of view of approximations. The problem 
is that the R.H.S. (right hand side) of (3.18) has momenta Pi and q effectively 
restricted to the range A — e < Pi, q < A + s (c.f. (3.14) and discussion below (3.4)), 
but either from (3.18) directly or by thinking about the correspondence (3.13) we see 
that the vertices of Sa sensitively depend on the form of 9 e there. (It is also helpful to 
think about this perturbatively i.e. in terms of connected Feynman diagrams using 
(3.11): the 6> e 's appear in every propagator). Of course the exact <S(pi, • • • , p n ; A) 
is not sensitive to the form of 9 e for momenta pi well away from A (and Ao), but 
it should be clear that if we attempt to approximate equations (3.18) directly, by 
truncation and/or otherwise, this property of the equations will ensure in general 
that our 'approximate' solutions depend very sensitively on the choice of 9 e and 
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truncation etc. The problem becomes especially clear if we consider the limit (3.5). 
Then the question is how to take the limit of the R.H.S. of (3.18) remembering that 
S(— Pi, h; A) for example is an as-yet unknown function of O e (Pi, A). The answer 
is not at all the usual physicist's expedient of putting 6(0) = \. 

To demonstrate this, let us prove a little lemma: 

<y e (p,A)/(0 e (p,A),A)-<y(A-p) fdtffap) as e^O (3.19) 

Jo 

where f(6 e ,A) is any function whose dependence on the second argument (A) 
remains continuous at A = p in the limit e = 0. This follows from the identity 

S £ (P, A)/(0 e (p, A), A) = | A jf^ A dt f(t, A') 

(c.f. (3.15)), by noting that (by the properties below (3.4)) the integral is a repre- 
sentation of a step- function in A but with height J^dt f(t, A'). As a consequence we 
have for example S £ (p, A)9 £ (p, A) — > |<5(A — p) as expected, but S £ (p, A)6 £ (p, A) — > 
J*(A-p). 

It follows then that we need to know precisely the dependence of the 5"s on 9 s (pi, A), 
where is any external momentum, in order to take the limit. On the other hand 
if we do not take the limit then we must be careful in approximations to accurately 
describe the dependence of the 5"s on the e (pi, A). For these reasons we need to 
resolve the limit by finding general solutions for ^(pi, • • • , p n ; A) parametrized in 
terms of functions that are continuous at pi = A. 

Actually, finding such general solutions is straightforward. For simplicities sake we 
specialize to the case that S\ has only even powers of (p, so that the Z2 symmetry 
<p <-> — <p ensures the S^pi, • • • , p n ; A) with odd n vanish. We also take A _1 (p) =p 2 
in (3.1), putting the bare mass term ^rriQip 2 in Sa [</?]• (Of course the results are 
independent of this split as we will shortly see: we choose it only because it makes 
the equations a little neater). 

Consider first the equation for the two-point function. From (3.18) this is 

^S(p- A) = K A (p)S 2 (p; A)-\ J AA K A (q)S(q, -q, p, -p; A) 
with b.c. S(p; A ) = m\ , (3.20) 
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where b.c. means boundary condition. This has solutions S that, while discon- 
tinuous at A = p in the limit e = 0, neither vanish nor diverge there. For such 
solutions 

where finite means a term that is finite at A = p in the limit e = 0. Integrating, 
using (3. 14) (3. 15), we thus have 

^ = £[».fcA)-».feAo)] + ^ (3.21) 

where E(p; A) is continuous and non- vanishing at A = p, in the limit e = 0, and 
satisfies the b.c. S(p; A ) = TOq. S(p; A) has the interpretation of an effective self 
energy, indeed in the limit e = we have 



S(p; A) = E(p; A) for p < A, 

L 

(3.22) 



while from (3.13) we have A) = — — -!- — — for A > p > A 

p 2 + T,(p;A) 



The latter equation also implies that £ is 1PI (one particle irreducible). 
Now consider the equation for the 4-point function. From (3.18) this is 

d 4 
— 5(pi, p 2 , Ps, P4; A) = 5(pi, p 2 , p 3 , P4; A) ^ K A (p l )S(p l ; A) 
aA i=i 

1 /■ 



1 f do 

2 y 727w KA ^^^ q ' -q ' pi ' p2 ' ps ' p4; a ^ 



(27T) 

thus 



d 4 

— In 5(pi, p 2 , P3, P4; A) = K A (pi)S(pi] A) + /inite 
aA i=i 

4 a 

= X! ^ n *^Cp*5 A) + finite . 

where in the last line we used (3.20), and finite now means a term that is finite at 
A = pi i = 1, • • • , 4 in the limit e = 0. Integrating up we have 

S(pi, p 2 , Ps, p 4 ; A) = T( Pl , p 2 , p 3 , P4 ; A) J[ ^ ^ (3.23) 



i 1 ! Sfe; A) 
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for some 4-point function r(pi, p 2 , P3, P4; A) which is continuous, in the limit e = 0, 
at the points A=pi i = 1, 2, 3, 4. The product over S's is allowed because it also has 
this property. It was slipped in because it gives V a neat interpretation, namely it is 
the 4-point vertex of both the Wilsonian effective action and the Legendre effective 
action (and thus 1PI). Indeed in the limit e = we have, using (3. 13), (3. 21) and 
(3.22): 

S(Pi,P2,P3,P4; A) = r(p 1 ,p 2 ,p 3 ,p 4 ; A) all Pi<A 

4 

C(pi, P2, P3, P4; A) = -r(pi, p 2 , P3, P4; A) Y[ G(Pi', A) all Pi>A . 

At this stage we take the hint that, by judicious parameterisation, the continuous (in 
fact smooth) parts of S^pi, • • • , p n ; A) can be identified with the 1PI (one particle 
irreducible) vertices of the Legendre effective action. From the structure of the flow 
eqns (3.18), we see that S(pi, • • • , p n ; A) has an expansion in tree diagrams which 
on any 1PR (one particle reducible) leg, carrying momentum J2 Pi ei Pi where I is a 
subset of the momenta {pi, • • • , p n }, there is a term containing 6> e (| J2 Pi ei p* | , A): 
it is this that gives the first term in (3.18). Indeed it is obvious diagrammatically, 
through the identification below (3.11), that «S(pi, • • • , p n ; A) is made up of 1PI 
bits which are continuous in the limit (since all 6> e 's are integrated over), connected 
by 1PR legs through terms containing the factor A>(| J2 Pi eiPi\)- We thus have 
immediately that, in the limit e = 0, only purely 1PI contributions remain if A > 
an I T, Pi eiPi\ f rom (3-13): 

S(pi, ■ ■ ■ ,p n ; A) = T(pi, • • • ,p n ; A) if A>|^p,| for all /C {pi,---,p n }, 

PiEl 

(3.24) 

where T(pi, • • • , p n ; A) are the 1PI vertices of the Legendre effective action for the 
theory with I.R. cutoff A. (Compare (3.22) and the previous eqn.). 

Now let us give a full non-perturbative proof of these assertions. First we note 
that by (3.9), W\[ip < , J] = IiiZ\[lp < , J] is in fact the generator of connected greens 
functions even for <p< ^0. Thus Ta defined by 

i(^ c - y^.A; 1 .^ - p<) + r A [^ c ] = -W A [¥><, J] + J.^ C (3-25) 
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is, by the usual analysis, the generator of 1PI greens functions. In here cp c is the 
classical field (p c = 8Wa/SJ. The extra terms on the left hand side are present in the 
classical action (c.f. (3.9)) and would normally be included in the Legendre effective 
action but are not strictly speaking 1PI diagrams, therefore we have written them 
separately. The reader may wonder if we are being sloppy here by not displaying the 
dependence of Ta on <p < . In fact it is clear that T\ has no perturbative dependence 
on cp < , since using (3.9) one cannot draw 1PI Feynman diagrams connecting cp < . 
Non-perturbatively it is true too since we have from the definition (3.25): 



5W A 



5(f < 



5p< 

but the R.H.S. vanishes, since from (3.11) one proves 

<v< > v « > 

Substituting (3.11) into (3.25), and rearranging to get $ = A>. J + </?<, one finds 
the following generalised Legendre transform relation between S A and T A : 

S A [$] = T A [ip c ] + i(^ c - <f>).A-\(<p c - $) (3.26) 

which in particular implies (p c = $ — A > .(5S A /5&). Substituting this back into 
(3.26) gives: 

By iteration, starting with S A [$>] = Ta[$], one obtains the now expected expansion 
of S A in tree diagrams (as follows from (3.11) and described above (3.24)). On the 
other hand for the vertex S(pi, ■ ■ ■ , p n ; A), if we ensure that A > | J2 Pi ei Pi I f° r 
all subsets I C {pi, • • • , p n } and take the limit e — > 0, then the above A> terms 
all vanish, proving (3.24). Physically, recalling that tree diagrams correspond to 
a small field expansion of classical field theory, we may interpret result (3.27) as 
follows: Ta is the Wilsonian quantum effective action obtained from integrating 
out purely quantum field excitations with invariant momentum p > A, while in S A 
in addition we integrate out also classical field excitations with momentum p > A. 
For sufficiently small external momenta a given tree level process with intermediate 
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momentum p > A is not possible, so here the vertices of the two effective actions 
coincide. 

We now derive the flow equations for Fa, from which we will see explicitly that 
the dependence on A is smooth, and the limit e = may be taken unambiguously. 
Substituting (3.11) into (3.16) we have 



dA 



J 



-V - V?<)-^V - <P<) - -tr ( 
2 1 ^ V<} dA { ^ 2 \ dA 5JSJJ 



which, using (3.25) and the relation 

S 2 W A = / ! 5 2 T A 
5J5J V > S<p c 6<p c 

(derived in the standard way from (3.25)), gives 



We now drop the uninteresting field independent part (i.e. the vacuum energy) from 
both sides; we must in any case for consistency since we have done so before (c.f. 
(3.2) and below). This requires separating from the above two-point function the 
(field independent) effective self-energy: 

; "' A X xy + t xy [<p c ] . (3.29) 



Hence, using the fact that E = S(p; A) is diagonal in momentum space and (3.14), 
the subtracted version reads: 



d_ 
8A 



^M = 4 tr {(TT^- f -( 1 + ' A > 1 + s '"- f )" 1 } 



(3.30) 



This is the version of the 1PI flow equations with a smooth cutoff. We see that 
there is no longer a term with an unintegrated K\, like the first term in (3.18), 
responsible for the discontinuous behaviour of S^pi, • • • , p n ; A) in the sharp cutoff 
limit. Indeed, since all 6>-f unctions now appear under the integral (trace), and they 
appear in such a way as to have a sensible limit, the R.H.S. is a smooth function 
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of A in the limit e = 0. The boundary conditions at A = A follow from (3.4) and 
(3.27): Ta = Sa = bare action. 

We may take the limit, using our little lemma (3.19) and the explicit relation (3.4) 
and (3.15). We specialize to A _1 (p) = p 2 as before (c.f. above (3.20)). This gives: 



8A 
where now 



d v A ^} = 1 fd ° q 6{q - A) 



2 J (2n) D q 2 + Z(q;A) 

0(p-A)-0(p-A o 



r^l + G.r)" 1 (q,-q) 



G(p; A) = 



(3.31) 



(3.32) 



p 2 + Z(p;A) 

The assumption made in using (3.19) is that the term in (3.31) enclosed in 
square brackets, or equivalently r[</? c ](q, p), is not singular — is in fact continuous— 
at p = — q. This is true providing ip does not have a non-zero constant part 
<f> (i.e. vacuum expectation value, in momentum space a term of the form 
(fi(p) ~ (27r) r> 5(p) <(p>). In the latter case we must split off from r[</? c ], the 
constant field part i.e. instead of (3.29) above one writes: 



s 2 r. 



8ip x 8ip c y 



(«p>) + r xy [<p c 



(3.33) 



where T, xy (<</?>) is still diagonal in momentum space and defined by 



H X y(<ip>) 



5 2 T 



A 



(p c = <<fi> 



8ip x 8ip 

Working again from (3.28) and dropping a field independent term we obtain 

dT A [<p' 



(3.34) 



<9A 



1 f dD( l s , ^ f 1 



r^l + G.r)" 1 ] (q, -q) +Vln(q 2 + E)j . 

(3.35) 



where now £ is given by (3.34) i.e. 
£ = £(<?; <y?>, A) and 



G(p; A) 



0(p-A)-0(p-A o ) 



p 2 + X(p;«p>,A) ' 
and V is the space-time volume. 

Using (1 + G.f) -1 = 1 - G.f + (G.f) 2 in (3.31) and expanding in tp c one 

obtains the flow equations for the vertices: 



d f d D 

M r(pi '---' p - A)= y?2v 



d u q 8(q - A) 



(2n) D q 2 + Z(q;A) 
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£(q,pi, • • • ,Pn; A) 



(3.36) 



where 

E(q, Pl, • • • , Pn, A) = -^ F (q 5 -Q» Pl» • • • > Pn! A) 

+ ]T r(q,-q-P 1 ,/ 1 ;A)G(|q + P 1 |;A)r(q-P 2 ,-q,/ 2 ;A) 

- ]T r(q,-q-P 1 ,/ 1 ;A)G(|q+P 1 |;A)x (3.37) 
{h,I 2 },h 

r(q + P l5 -q + P 2 , 7 3 ; A)G(|q - P 2 |; A)r(q - P 2 , -q, 7 2 ; A) 

+ • • • . 



Here Pj = J2 Pk eii Pk and Z){/ 1 ,/ 2 },j 3 ,...,7 m means a sum over all disjoint subsets 
Ii Hlj = (Vi, j) such that U™i 7j = {pi, • • • , Pn}- The symmetrization {7i, 7 2 } 
means this pair is counted only once i.e. {ii, 7 2 } = {/ 2 , 7i}, and is a reflection of 
the Bose symmetry in the equation. Once again the equation is best appreciated 
graphically — c.f. fig. 2. 

Evidently the expansion stops at the term where all 1PI vertices have their minimum 
number of legs, i.e. at the n th term in general, or the (n/2) th term in the (p <-> —(p 
invariant theory. It is often helpful to rewrite (3.36) by performing the integral over 
q and factoring out the 7>-dimensional solid angle thus: 

Ar( Pl , • • • , p n ; A) = T{ ^ {Aii) o A 2 + E(A;A) < Pi, • • • , p.; A) > g . A . 

(3.38) 

In here < • • • > q= \ means an average over all directions for q, its length restricted 
to q = A. 

Of course these equations can be derived by first expanding the smooth cutoff 
equations (3.30) and then taking the limit. Equations (3.36) (3.37) are valid 
providing the external momenta in the propagators G do not vanish {i.e. Pi, P 2 7^ 
etc. ). If this happens then, by (3.32) and (3.36), we obtain 0(0) 's which, as we 
have already seen, are ambiguous. If nothing special happens at zero momentum 
we can ignore this problem since they are points of measure zero in any calculation 
or realistic physics problem. On the other hand zero momentum is special if the 
field acquires a vacuum expectation value; the correct equations either follow from 
(3.35), or expansion of (3.30) and the setting of the relevant momenta to zero before 
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Fig.2. The flow equations for the 1PI vertices in the sharp cutoff limit. Internal 
lines are full propagators. The black dot now represents restriction to momentum 
q = A; the other propagators have an I.R. momentum cutoff p > A. 

using the limit (3.19). Apart from the odd aside, we will use (3.36) (3.37) from now 
on and ignore the more general case (3.35) in this paper. 

Equation (3.32) implies that a contribution in (3.37) vanishes if an intermediate 
momentum, e.g. |q + Pi|, is larger than Ao and so has a profound effect on 
the solutions especially for external momenta and/or A close to (but less than) 
the cutoff. If one recalls however that T\ may be regarded as just a smooth 
parametrization of the effective action S\ this result might seem a little puzzling: 
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after all in the Polchinski equation (3.17) contributions from momenta p > A + e are 
heavily suppressed, and in the limit e = one naively would expect the equations 
to depend only on momenta p < A. As we have seen however, although all the 
£(Pi> • • • > Pn! A)'s are restricted to pi < A in the limit e = 0, in the sense that any 
propagator connecting them is A K (p) = for p > A (c.f. (3.7) and (3.10)), the limit 
itself is ambiguous. To resolve the limit as the cutoff becomes sharp it is necessary 
to isolate the rapid change at the cutoff from the smooth parts. The latter enter the 
limit in such a way that they need to be known at least as a series expansion (about 
Pi = A) to all orders. In other words we need a continuation of the smooth parts to 
the region pi > A. There are presumably many ways to extend the smooth parts to 
Pi > A but nota bene however that the resulting S(pi, • • • , p n ; A)'s, in the regime 
Pi < A, will be invariant under such choices. If we require that these choices also 
correctly represent 1PI greens functions for pi > A then the number of possibilities is 
severely restricted. The simplest possibility is to require that the r(pi, ■ ■ ■ , p n ; A)'s 
are taken to be as analytic as possible, and the minimal requirement of the series 
expansion about Pi = A incorporated as the maximally analytic extension into the 
region Re(pi) > A. This is essentially what will be assumed in our approximation 
scheme. It is clear that such a prescription involves dropping the A terms in (3.4) 
and (3.32), in other words taking the limit A — > oo, since any Taylor expansion 
about p = A will not see the theta function 6{p — A ). 

One can see that no U.V. singularities are introduced into the eqns.(3.30)-(3.38) 
by the limit Aq — > oo since they involve only an integration over q with fixed 
length q = A. Of course U.V. divergences will still arise in the solutions, so we 
need now to set the boundary conditions at some finite scale A = Aq. The most 
analytic choice will be obtained if the boundary values r(pi, • • • , p n ; A ) are chosen 
to be analytic functions of momenta , for localities sake we further choose the 
'bare' [(p c ] to have a local derivative expansion. The relation between the 'bare' 
Wilsonian effective action Sa> and this 'bare Wilsonian quantum effective action' 
now involves the full tree expansion (3.27). One can expect though, by the usual 
universality arguments, that these differences only alter the way the cutoff effects 
are parametrized in any given theory (compensatable by a reparametrization of the 
irrelevant couplings), while the universality classes remain the same. However, the 
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concept of a bare action in a path integral as in (3.1) appears now to be lost. 

Actually nothing quite so drastic has happened: we know that we can reconstruct 
the theory for p < A' by using S\' as the bare action and an overall momentum 

A' 

cutoff Aq, and thus also a Legendre effective action r A ° if we use e.g. (3.31) and 

A' 

(3.32) with Ao i— > Aq and the boundary condition T A ? = Sa' - We display explicitly 
with a superscript the fact that this new Legendre effective action has an overall 
cutoff Aq. From the boundary condition if nothing else, it evidently differs from 
the r A ° — > discussed in the previous paragraph. By construction however it 
nevertheless yields the same S\ for all scales less than A . We have thus sketched 
how to recover a Legendre effective action computed from a bare action at A by 

A A' 

a reparametrization of the Legendre effective action r A ° i— ► r A ° , while leaving the 
physics at scales lower than A untouched. This is thus an invariance of the form 
stressed above. 

We now fill in the details. As already indicated in the introduction the reparamet- 
rization is itself just a tree diagram expansion. Subtracting the field independent 
term ±tr{(dA>7 rfA )- A >} from (3.28) we have that 



<9A 



1 

-tr< 



dA 



§ 2 T Ao 



A 



8ip8ip 



(3.39) 



In here we have explicitly displayed the U.V. and I.R. cutoffs on A> 



with 



Aj° = [0 e (p,A)-0 e (p,Ao)] A(p) 



(3.40) 



as given in (3.4). We have dropped the superscript c on the classical fields to 
unclutter the notation. Define the required tree diagram relations through the 
generalised Legendre transform: 



(3.41) 



A A ? is defined as in (3.40), and here as elsewhere A < A . The tree expansion for 

A' 

r A ° follows by solving for ip iteratively via 



A I 

5ip 



+ (Ajj)" 1 .(^- v /) = , 



(3.42) 
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as 



(p = (p 



-i <v 

' 6<p' 



+ 



(3.43) 



and substituting back into the R.H.S. of (3.41). Here E A ° is the self-energy derived 
from the 2 nd derivative of rj°. (Note the close similarity of (3.41) and (3.26). The 
analogous expansion to (3.43) gives the relation between S\ and F\ in terms of full 
propagators.) 

It is instructive to perform this for low orders and see diagrammatically how the 
tree expansion simply shifts the overall cutoff in the equations (3.36) (3.37) (3.32). 
Here we will confirm it in full, analytically. From (3.41) we also have 



Sip 



fcl + (Ajj)- 1 .(„-*/)=0 , 



(3.44) 



which differentiating again with respect to <p' gives 



2r A () 



S 2 T 



Sip' I Sip'Sip' 



A 



yz 



Sip'Sip' 



- (Aj?) 



On the other hand, comparing (3.42) and (3.44) we have 5T%° /6<p' = ST%°/8(p which 
on differentiating again with respect to <p' gives 



S<p' \ Sip'Sip' 



\ S( P S( f J. 



Combining the last two displayed equations, and using A A ° = A A ? + A A ° as follows 
from (3.40), we have 



+ Aa° = 



s 2 rl° 

Sip'Sip' 



+ Aa° 



(3.45) 



A A' 

Together with the identity dA A ° /dA = dA A ° /dA, which trivially follows from (3.40), 
and the identity 



d_ 
OA 



r A° Vp'\ 
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which follows in the usual way from (3.41) and (3.44), we see that (3.45) implies 
that the flow equations (3.39) are unchanged by the reparametrization (3.41) i.e. if 

A A' 

r A ° satisfies the flow equations then so does r A ° (and vice versa). It only remains 
to show that physics with momentum p < A' is unchanged. But this follows because 
the new flow equations, (3.39) with Ao i— > A , have the correct boundary condition 

A' 

r A ? \(f'\ = Sa' q [<£>']• To confirm this last point set (p c = cp, <E> = <// and A = Aq in (3.26); 
recall that in (3.26) everything has overall cutoff A and compare with (3.41). The 
result follows immediately. 

4. Approximation Methods. 

In this section we consider methods of approximating the sharp cutoff flow equations 
(3.36)-(3.38) where we drop the overall momentum cutoff as discussed at the end 
of the previous section. We start however by making some general comments on 
the non-local effects induced by taking the sharp cutoff limit. We then turn to the 
flow equations with smooth cutoff: these do not suffer from non-localities, but have 
the serious calculational problems already mentioned in the introduction. After 
considering the straightforward numerical approach, we look at the most promising 
method, an expansion in external momenta about pi = in the sharp cutoff limit. 
This is briefly compared to two alternative expansions. 

The non-localities correspond to non-analytic dependence on the momenta at = 0. 
We have already seen an example of this in the analysis of (3.33)-(3.35) where the 
constant part of (p had to be treated separately. This is connected to the fact that 
the solutions to (3.36), equivalently (3.38), generally do not have well defined limits 
as any pi — > 0, rather they remain functions of the relative angles between the 
Pi even in the limit. One can see this by inspection of the third term of (3.37), 
however first let us look at the second term. In the greens function G(|q + Pi|; A) 
there is a term 0(|q+ Pi| 2 - A 2 ) ss 6>(q.Pi) for P x « A. Expanding the T's in 
the second term of (3.37) in this limit we see that there are terms of the form 
< q.Pi#(q.Pi) > 9= a~ A|F*i | , so that r(pi, • • • , p n ; A) has an expansion in lengths 
|p| and not p 2 as would be required by locality /analyticity. In the third term one 
also gets contributions of the form < 0(Pi.q)0(P 2 .q) > g =A when Pi, P 2 « A which 
depend only on the angle between Pi and P 2 and not on their lengths at all. 
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Disturbing as these non- localities might at first appear they are all understandable 
in the intuitive sense already covered in the introduction. Indeed one can see all 
these effects already in perturbation theory i.e. in Feynman diagrams, if one provides 
each propagator with a sharp I.R. cutoff (as in (3.32)). We have searched for any 
deep problems with using a sharp cutoff, this is a primary reason for analysing the 
simple models in the next section, but we have not uncovered any. One obvious 
constraint is that the resulting physics in the limit A — > does not suffer from these 
non-localities. In other words that the non-localities present in the intermediate 
effective action Ta are cancelled by those in the remaining physics with momenta 
constrained to have p < A. But this is guaranteed if we keep the bare Fa [</j] 
local. 8 (Presumably it is possible to be more subtle, noting that since the typical 
effective action Ta is non-local the bare effective action can be too, providing the 
non- localities have some closely prescribed form). It seems possible that non-local 
terms in Ta could turn out to be relevant — in other words have to be added as 
counterterms to the bare lagrangian, making it non-local, in order to renormalize 
the low energy theory. It is clear that this does not happen in perturbation theory 
from the interpretation in terms of Feynman diagrams above; it is also clear that it 
cannot happen non-perturbatively in any truncation, by thinking about the smooth 
case first and then taking the limit e — > (because the divergent parts of the 
bare couplings will not depend on s in this limit). But it is allowed by naive 
power counting. If this would happen as the result of some approximation in the 
momentum dependence (c.f. the expansions later), it would be a disaster: the low 
energy physics would likely be non-local also, while new and unphysical parameters 
would have been added to the bare lagrangian. As we will see in the next section 
this doesn't happen, at least in the situations we considered. 

It is intuitively clear that taking the sharp cutoff limit ought to simplify matters, but 
given the non-local behaviour sketched above it is worthwhile taking a careful look 
at approximating the smooth cutoff equations (3.30). This can be done by making 
a momentum expansion around pi = and truncating at some order. Recall the 
general form of £ as described below (3.4) . From there it is clear we must have e < A, 
otherwise there will be no effective cutoff on low momenta. Similar 9 terms as in the 

We now drop the primes introduced at the end of the last section. 
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sharp cutoff case are of course encountered but now they appear as e (|q + p|, A) 
where p = e.g. Pi, P2 etc. Expanding these about p = gives 

e (|q + p|, A) =6 £ (q, A) + ^f^- q 0e{q, A) 

1 [ pV-(g.p) 2 a A ^ (q-p) 2 <9 2 M \ 

(4.1) 

We see that because 9 e is a smooth function a true Taylor expansion in the about 
Pi =0 now exists for the r(p 1; • • • , p n ; A), equivalently T\[(f] has a local derivative 
expansion, with radius of convergence p e. This radius follows from (4.1) since we 
have by simple dimensional arguments (d/dq) n 6 e (q. A) p=* |e~ n for q ps A, and since 
e < A. Substituting the expansion of Ta into (3.30) however involves integrating 
the expansion over a momentum shell with radius A, thus for convergence of the 
expansion we need A < e. We see we need both e < A and e > A. We can relax each 
of these to equality but then we have only a weak suppression of low momenta and 
an expansion at its radius of convergence where it will converge only very weakly if 
at all. 

To justify the estimated radius of convergence we look at a few examples of 9 e . If 
9 £ is an exact cutoff on either side i.e. 6 e {q, A) = for all q > A + e and/or for all 
q < A — e then the radius of convergence of the expansion (4.1) with q ps A is at 
most e from elementary complex analysis. An example of an approximate cutoff is 

(9 e (g,A) = i(l + tanh[|(g-A)/e]) . 

A plot of tanh(3x/2e) will convince the reader that this has width ps 2e as required 
(tanh(3/2) = 0.9). On the other hand the radius of convergence of (4.1) is p = 
7re/3 ~ e as given by the first pole of tanh(^) on the imaginary z-axis. Actually it is 
technically possible to have infinite radius of convergence, but not in any practical 
sense since for e < A the series only starts to converge at high order. For example 
write 

9 £ (q, A) = \ (l + erf [a(q 2 - A 2 )/2Ae]) , (4.2) 

where erf is the error function and a is to be chosen so that (4.2) has width 2e. 
Clearly the convergence is best if we take the maximum width e = A. We require 
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a ss 2.5 to get effective suppression of low momenta (6>a(0, A) = 0.04), i.e. to get 
the correct width in the q < A regime. But the fact that q appears squared in (4.2) 
makes the q > A behaviour very sharp, and this has serious effects on the effective 
convergence of the series. Indeed estimating the convergence in the flow equations 
by expanding (4.2) about A as q = A + Aq, and then setting Aq = ±A, one finds 
that 10% accuracy requires expansion to 84 th order! 

The second major problem with the smooth cutoff case is as follows. Substituting 
(4.1) (and an expansion of F\) into (3.30), one sees that the coefficients of the 
momentum expansion for Ta are determined by integrals over q of 8 S (q, A) times 
terms containing (d/dq) n 9 s (q, A). Unlike the terms in (3.19), these coefficients 
have no universal behaviour for s « A, but depend sensitively on the shape of 
the cutoff. 9 While this is reasonable, given that these higher derivative corrections 
are 'mixing' with the 9 e (q, A) -1 in the kinetic term — the latter corresponding to a 
higher derivative regularisation with an infinite number of derivatives, the lack of 
universality makes their interpretation difficult. Furthermore, if we take e ~ A the 
sensitivity to the form of the cutoff is physical since it reflects how much of the low 
momentum modes q < A are integrated out. 

There are a few less major problems with the smooth cutoff case also: the above 
integrals can be cast as one-dimensional integrals over q but they still have to be 
done numerically, and since the integrals depend on E(0; A), (3.30) becomes a set 
of integro-differential equations. Also the imposition of a wide smooth cutoff makes 
matching to other methods — e.g. numerical, perturbative etc. — problematical. 
Given all these problems, it seems clear that local approximations to the smooth 
flow equations can not be extended in any substantive way beyond the uncontrolled 
approximation for the effective potential (covered in the introduction) where it 
differs inessentially from the sharp cutoff limit. 

Let us note here that attempts to solve either the smooth or sharp cutoff flow 
equations directly by numerical means will also encounter severe problems. We 
would need to discretize over the arguments of the greens functions, replacing d/dA 
by a finite difference and the integration over spherical shells by the appropriate 

Of course they diverge as e — > 0, but we mean also that the coefficient of the 
divergence has no limit. 
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sums. At least naively, we must work with the values of the r(pi, • • • , p n ; A) over all 
combinations of discretized values for the pf i = 1, ■ ■ ■ ,n fi = 1, ■ • ■ ,4 in the range 
< pi < Aq. Even allowing for rotational invariance, permutation symmetry, and 
momentum conservation, this quickly gets out of hand as n increases. The problem 
is not so much the storage of so many numbers but the fact that they all have to 
be updated as A i— > A — 8 A. It is easy to convince oneself that the computational 
difficulty is greater than that of lattice gauge theory even if we consider only n < 6. 
There seem to be many possibilities for improving this naive scenario but the crucial 
observation is that much of the integration domain in (pi, • • • , p n ; A) space must be 
boring: this suggests that we should try to handle these regions semi-analytically, 
e.g. by the momentum expansion which follows, leaving numerical analysis for 
(hopefully) small ranges where semi-analytic methods may fail. Actually, semi- 
analytic methods could be pushed very far before their numerical complexity rivalled 
that of straightforward numerical approaches. 

We turn at last to the most promising approach. We work with the sharp cutoff 
equations (3.37) (3.38). Writing all mass scales in units of A utilises the scaling 
symmetry in the equations. Thus we replace A with t = ln(A /A) (initial b.c's 
are set at t = 0), and write T(Api, • • • , Ap n ; A) = A d ™7(pi, • • • , p n ; t) where d n = 
D + (1 — D/2)n is the engineering (or canonical) dimension of T(pi, ■ ■ ■ , p n ; A). 
Similarly we write E(Ap; A) = A 2 cr(p; t). Finally we replace pi by ppi where p is 
the 'momentum scale' and is our expansion parameter. (It may be set to one at the 
end). This gives 

( d d , \ , 

I -Q t + Pq- - d n pyppi, • • • , PPn; t) = 

2 1 (4 ' 3) 

' (4,)fr(f)l + a(l;i) <% ^-'^ f)> ^ • 

E is still given by formula (3.37), (with T's replaced by 7's) see also fig. 2, while G 
is now given by 

n( n - +\ = 

p 2 + a(p;t) 

Now we expand in small momentum scale as follows 



*)-^ ■ (4-4) 



7(PPl» • • • » PPn] t)= P m 7m(Pl, ■■■,Pn,t) (4.5) 

m=0 
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so the 7 m scale homogeneously: 

7m(pPl,---,PPn;t)=p m 7m(pi,-",Pn;t) • (4.6) 

This expansion incorporates the non-local effects covered at the beginning of this 
section since the integer m is not restricted to being even and e.g. 70 is not 
necessarily constant but is in general a function of the angles between the p^. The 
R.H.S. of (4.3) may be evaluated as follows. Let x = p.q/p = cos $, where d is the 
angle between q and p, then expansion of the I.R. cutoff is achieved as 

9(\pp + q| - 1) = 9(x + pp/2) = 9{x) +E ^| (| )" S^~ l \x) . (4.7) 

Terms of the form 7 m (q + pPi, — q + /0P2, • • •) (the most general) have re-expansions 
in p. The angular average may be performed, for example in four dimensions, 
writing q = (cos sin# n) where n is a unit 3- vector: 

2 /-7T 2 f 1 

< • • • >o=i= - / dd sin 2 # < ■ ■ ■ > n = - dxVl - x 2 <•••>„ . (4.8) 

7T JO - 7T J-l - 

The R.H.S. has an average over azimuthal directions. The approximation is achieved 
by working to some maximum order in p for each n-point function, and to some 
maximum n e.g. by setting higher n-point functions arbitrarily to zero. Since (4.3) is 
a set of first-order differential equations with boundary conditions fully determined 
by the 7(pi, • • • , p n ; 0) of the 'bare' action, the solutions to these truncations are 
unique. In the cases we discuss in the next section the solutions for 7 m factor into 
simple functions of the momenta times t dependent coefficients. 

This approximation can be expected to work well if rapid variations with respect 
to momenta happen not at all, or at a scale p such that p/A « 1. In addition if 
the 7's scale approximately as ~ A~ d ™ (i.e. T's ~ constant) below some mass scale 
M then the expansion is really in p/M. Mathematically one expects the radius 
of convergence to be the maximum one determined by the first singularity in the 
complex p plane. 

As discussed already in the introduction the crucial simplification in this approxi- 
mation scheme is the expansion of the I.R. cutoff (4.7). There are a couple of other 
realistic systematic expansions, but they are probably of limited use. One is to 
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reintroduce an overall cutoff A : After scaling out A, everything may be expanded 
in t = ln(A /A). Thus the cutoff terms (c.f. (3.32)) become: 

B(p - 1) - 0(p - e *) = tS(p - 1) + | [S(p - 1) - 5\p - 1)] + • • • (4.9) 

and 7(pi, • • • , p n ; t) = J2m=o ^ m 7m(pi 5 • • • , Pn)- The r 's are the bare vertices. 
Evidently this scheme only works for A not much different from Aq. If Aq is the 
true cutoff such an expansion will only see cutoff effects, so it is only useful if Aq is 
an effective cutoff, e.g. introduced via the reparametrization discussed at the end of 
sect. 3. The expansion (4.9) introduces via the angular average (where p is replaced 
by |Pfc + q|) powers of 1/Pfc. This expansion does not therefore converge for small 
momenta and does not match easily on to the above momentum expansion. 

Another possibility is (again) to drop A , not scaling out A but working with (3.37) 
and (3.38), and expanding in A « p t . This implies also Taylor expanding in q^. In 
this regime we can set all terms to 1. The angular average is easily computed by 
the identities < q^q u > q= A= 5^A 2 /D and their generalisations. This expansion is 
easy to perform and works well in the, albeit limited, regime A « pi. 

5. Examples. 

In this section we investigate two simple models of four dimensional \(p 4 theory to 
test the approximation scheme (4.3)-(4.8). The first is the analogue of the ladder (or 
cactus) approximation used in Dyson- Schwinger equations: it allows a comparison 
and checks on the effect of a sharp cutoff. The second model incorporates the 
first non-perturbative correction to the coupling and provides some tests of the 
momentum expansion. We also check that no non-local terms become relevant. An 
overview of the results has already been given in the introduction. 

To obtain the cactus approximation to the flow equations (4.3) we truncate to the 
n = 2 equation: 

di + P dp ~ 2 j ^ t] = W 1 + a(l; t) < 7(q ' " q - P ' " P; t] ^ • (5J) 

Now we substitute the ansatz 7(pi, • ■ ■ , P4; t) = A which corresponds in particular 
to ignoring the running of the coupling. It also implies in this approximation that 
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a = cr(; t) has no momentum dependence. Thus, writing a = X/ (4tv) 2 : 

While we ignore the running of a in this equation, we keep in mind that as we flow 
to low energies (i.e. increasing t) a shrinks. 10 The shrinking will continue to a = 
at A = (i.e. t = oo), the gaussian fixed point, unless a becomes bigger than one 
below some scale A = M. In the latter case we expect that a freezes out at a value 
~ a(M) for scales A < M (while £(; 0) ~ M 2 ). Keeping in mind this qualitative 
behaviour will ensure that we interpret the results of (5.2) correctly. 

Consider first the case where the physical mass-squared £(; 0) > 0. In this case the 
solution to (5.2) is easily figured out: 

(5.3) 

where ao = a(; 0) = ttiq/Aq. 

Equation (5.2) can be recast as an integral equation in a suggestive way. It is easiest 
to go back and integrate both sides of (3.36) between A and A = Ao, again with 
n = 2 and T(pi, ■ ■ • , P4; A) = A: 

A r A « d D q 1 
I a (2^p ^£(; q) 

The range of the integration is A < q < A as indicated. For A = this is the cactus 
approximation to the Dyson-Schwinger equations, except for one alteration: the 
self-energy inside the integral is not £(; 0) but is evaluated at scale A = q. One 
might expect that such an alteration gives a better approximation. The cactus 
approximation to the Dyson-Schwinger equations gives a constant £, determined 
implicitly by: 

(5.5) 



/ . s 9 A f A ° d D q 1 . , 



2 a 

m ° + 2 



Ag-£lnfl + M' 



To get some idea by how much the two approximations (5.3) and (5.5) differ, we 
compare them at strong and weak coupling. At weak coupling they agree to order 



10 as follows e.g. from the 1-loop /3-function, but can also be seen qualitatively directly 
from the n = 4 equation. See later. 
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a (the result obtained by putting E = tjIq in the R.H.S. of (5.5)). This is simply 
because they both correctly incorporate the one-loop correction. More non-trivial 
is agreement up to a factor of V2 at strong coupling. For the flow equations: 

E fa Ti _ / 1 \ 
For the Dyson-Schwinger equations: 

Since the two cactus approximations agree at small a and differ by 40% as a — > oo 
they can be expected to be in rough agreement for all a > and for all reasonable 
values of ttIq/Aq. They differ significantly when one starts to try and tune E « Aq. 
This is covered later. 

To better appreciate the origin of the difference between the two approximations, 
iterate the equation (5.4) to get the perturbative expansion. The result is given in 
terms of Feynman graphs in fig. 3. 




Fig. 3. Diagrams contributing to the cactus approximation to the flow equations. 
The loops have the I.R. cutoffs shown in their centres. 
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It is the same graphical expansion as that of the cactus approximation in Dyson- 
Schwinger equations except that part of the momentum integrations are missing 
as a result of I.R. cutoffs. For A = the disagreement shows first at two loops. 
The remainder of this diagram is provided by adding the one-loop correction to the 
4-point vertex (in the ansatz above (5.4)) as indicated by the dotted box. (This 
gives the same diagram except that now the lower loop is I.R. cutoff by the higher 
loops momentum. Exchanging the order of integration here gives a contribution in 
which the upper loop is U.V. cutoff by the lower loops momentum, so the addition 
of the two contributions gives the two-loop diagram with momenta integrated over 
the whole range.) Of course if we incorporate the one- loop correction to the 4- 
point vertex then, by the symmetry of the flow equations we also generate the 
other two-loop self-energy graphs, as well as parts of many others. In contrast the 
Dyson-Schwinger equations organise the graphical expansion topologically; the term 
neglected in the Dyson-Schwinger equation for the self-energy generates all graphs 
in which the two external legs are joined to separate vertices. 

It is clear from (5.3) that the cactus approximation to the flow equation has a 
different behaviour for a < \ compared to a > \. The latter can be obtained 
from (5.3) by analytic continuation, noting the reality of the solution, but the 
qualitative nature of the solutions is best appreciated by studying the R.H.S. of the 
flow equation (5.2) as a function of a. This is sketched in fig. 4. 

In this figure a± are the zeroes of the R.H.S: a = a± = — \ ± \\J\ — 2a. For a < \ 
it is clear that if we require £(; 0) > we must have o"o > ot+. In this case as t — > oo 
we have a — > oo such that da/dt — > 2a i.e. E(; A) — > a positive finite constant, as 
A — > 0. The constant is of order Aq for generic <j$ (e.g. by dimensional arguments). 
Since a + is an I.R. unstable fixed point we can tune to £(; 0) « Aq by taking 
ml — > a + Ag; indeed if TOq = a+Ag then £(; 0) = 0. This behaviour is similar to 
that of the Dyson-Schwinger cactus approximation (5.5) where, for all a, S may be 
tuned to zero by ttiq — > — |oAq. 

For a > \ however, there are no finite fixed points in fig. 4. Now S(;0) ~ Aq 
or larger for all choices of m^. ( The phase given by < — Aq is covered 
below). This behaviour is not modelled by the cactus approximation to the Dyson- 
Schwinger equations. Is it reasonable? Keeping in mind the comments below 
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Fig. 4. da/dt as a function of a in the cactus approximation for (a) a < h, (b) 
a > \. Direction of flow for increasing t is indicated by arrows. 

(5.2) we see that we should interpret a = \ as an effective maximum renormalized 
coupling A max ~ 87r 2 = 79.0, above which there is no continuum limit because the 
renormalized mass is of order the cutoff. This is therefore already a reflection of 
triviality of four dimensional A<£> 4 theory. The accurate analysis of Liischer and 
Weisz[ll] gives A max = 78 ± 3 (!). (See table 3 of that paper. Note m\ = £(; 0)/Ag 
and gR = \). Of course the exact agreement is an accident since the cactus expansion 
is only a crude approximation, and uses momentum cutoff, not lattice cutoff as in 
ref.[ll]. 

Now consider the phases where cr(; t) flows to negative values. For a < \ we have, 
by fig. 4, an I.R. stable fixed point cr(; t) — > ct_ as t —>■ oo, for all bare masses in the 
range — 1 < <7q < a+. Thus the physical mass £(; 0) = for all bare masses in this 
range. If now we again recall the qualitative evolution of a given above (5.2) we find 
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reason to doubt this continuum limit. Indeed we expect that a — > as A — > in this 
regime, but this would cause the fixed point ct_ — > —1, the tachyonic singularity 
in fig. 4. From our general comments on convergence below (4.8), and explicitly 
from the factors (4.4) in the expansion (3.37), one can see that our expansion, of 
which the cactus approximation is the lowest order, would then break down. To 
understand what really happens for this range of bare parameters it is necessary to 
go beyond this simplest approximation. For example it is clear that the problem 
is avoided by first shifting the field ip e.g. to the minimum of the potential, but 
it may be that the problem is resolved after a reparametrization or as a result of 
unexpected behaviour in the flow of A. 

Finally, for all a > we have that Wq < — Aq yields a finite 'physical' mass 
£(; 0) < — Aq. Now however, expanding about the unstable symmetric point (p = is 
questionable independent of our approximation: Since the negative mass-squared is 
greater than the cutoff- squared, the symmetric phase cannot exist in any physical 
sense. In fact in this regime (5.2) is not the limit of the equations with smooth 
cutoff because the continuity condition in lemma (3.19) is violated for some q < A . 
If one figures out from (3.30) the smooth cutoff form of the cactus approximation 
then the R.H.S. is divergent if cr(; t) < —1, for all e > 0. 

The next step in the approximation method is to expand in momentum scale. We 
see from (5.1) that the n = 4 flow equation is needed to determine the coefficients. 
From (4.3) we obtain 

(I +P *W) 7(Pl ' ' ' • ' P4; * } = (4^*1 + *(!;*) < 7(q ' _q ' Pl ' • • • ' P4; t] >q ^ 
2 1 

~ tt^2 ~r, — rr^ < G ( r i2; *h(q» - r i2, Pi, P2; t)^{r 12 , -q, p 3 , p 4 ; t) 
{Any 1 + ct(1; t) 

+ (2 ^3) + (2 ^4) > q=1 

(5.6) 

where r ^ = q + Pi + Pj , and in the second average the last two terms are the same 
as the first term but with indices swopped as indicated. To truncate at this level 
we set 7(pi, • • • , P6; t) = 0. This is a sensible approximation because, even if the 
6-point vertex is non-zero in the bare action, it is irrelevant at the gaussian fixed 
point and therefore at sufficiently low energy-scales it shrinks rapidly as A falls, 
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until it is determined to good approximation as a perturbative series in the now 
small renormalized coupling A. Note that since the R.H.S. of (5.6) is positive in 
this truncation we can see already that the 4-point vertex must shrink as we flow 
to low energies. 

We now make the theory massless by setting a(p; t) = 0. This is not meant to be 
taken seriously as an approximation: we do so purely because it greatly simplifies 
the equations since we can now study (5.6) without having to worry explicitly about 
(5.1). With these changes we have 

< X \ 7(q> ~ r i2, Pi, P2; *h(ri2, -Q» P3, P4; t) (5.7) 



(4tt) v 12 

+ (2^3) + (2^4)> g= i . 

Determining the perturbative solution, we will find that this truncation is exact 
at one-loop. Again the perturbative expansion is easiest performed by integrating 
both sides of (3.36) with respect to A using the b.c. T(pi, • • • , P4; Ao) = Ao, and 
iterating. The Feynman graphs up to two-loops are given in fig. 5. 



Fig. 5. Feynman graphs contributing to the 4-point vertex, in the truncation defined 
by setting the 6-point vertex to zero (with the self-energy also set to zero). 

The one-loop contribution to T(pi, • • • , p 4 ; A) is found to be 

_ A 2 v fd 4 q [9(q - A) - 6(q - A )]g(|q + Pl + Pi \ - g) 
°hh^) 4 4 2 (q + Pi+P,) 2 

Substitute q — > q — p followed by q — > — q and add the result to the above (dividing 
the whole by 2). After a little rearrangement this is seen to be 

_ Ag * f d 4 q 9(q - A)9(\q + Pi + Pi I - A) - 9(q - Ap)g(|q + P i + pj - Ap) 
2 hh^Y 9 2 (q + Pi + P,) 2 
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This is the exact one-loop answer, with however a U.V. cutoff that only works when 
all momenta in the loop are larger than A . This is a reflection of our dropping the 
explicit overall cutoff at the end of sect. 3. Similarly the first two- loop diagram in 
fig. 5 is given exactly (as follows from adding several contributions after exchanging 
orders of integration) , while some of the second is missing because it is provided by 
the one-loop 6-point vertex (as indicated by the dotted box). 

Since eqn.(5.7) incorporates the exact one- loop result it gives in particular the right 
one-loop divergence — 3Agt/(47r) 2 . (Recall that t = ln(A /A)). If the truncation is 
indeed perturbatively renormalizable (c.f. sect. 2) then the leading log divergences 
~ XQ l t m ~ 1 , for all m, must be those determined by the one-loop renormalization 
group. Clearly these leading log terms come from iterating the constant part 
of r(pi, • • • , p4; A), i.e. dropping all external-momentum dependent pieces. But 
substituting a constant T(pi, • • • , p 4 ; A) = X(t) in such a way into (5.7) gives: 

i.e. the one-loop /3-function as required. 

It is an instructive exercise to compute the two-loop divergences. This will provide 
further explicit checks on the method. Two-loop diagrams of the first type in fig. 5 
are straightforward and give in total 3A[jt 2 / (4-7r) 4 . The last type of two- loop diagram 
gives in total 

rA d 4 p 1 r A o d 4 q 6>(|q + p|— q) 

OA, ' 



°Ja (2tt) 4 p 4 J p (2tt) 4 (? 2 (q + p) 2 ' 

(We have set pi, • • • , p4 = since they do not appear in the divergences). The inner 
integral is not doable exactly. We follow the momentum scale expansion (4.5)-(4.8) 
which here means expanding the inner integrand in p/q and performing the angular 
integration. One finds 

. / / \ 3 / \ 5 

A o d^q #(|q + p| — q) 1 f A °dq I tt lp 1 I p\ 3 I p\ 



r A o d 4 q ^(|q + p| — q) 1 f A °dq J 7r lp 
Jp (2tt) 4 Q 2 (q + p) 2 = 4^J P ~q~\4~6q 



240 \qj 8960 \q J 

(5.9) 

which, apart from the constant term, is an expansion in odd powers of p (as follows 
from taking the symmetric part under p — > — p and using < ^q^p > q =q> p = 
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Integrating, and substituting into the outer loop integral one obtains the O(Aq) 
contribution to r(pi, • • • , P4; A) as 

9A| / A \ 2 Ag A 

where 77 is given as a rapidly convergent numerical series 77 = ^(8 + ^ + +'"')• 
The partial sums 77 = 2.5465, 2.5677, 2.5687, • • • , 2.568818 converge to 3sf already at 
order p 3 in the momentum expansion, after which approximately an extra decimal 
place in accuracy is added with each new term. (The "finite" term converges 
somewhat faster). Not a bene that, while the terms bounded by curly brackets 
in (5.9) may be regarded as an expansion in the small quantity p/q, once the 
integral is performed there are no longer identifiable small parameters to control 
the approximation (because the lower limit is p) . Nevertheless, as we have seen, the 
expansion converges very rapidly. The leading log divergence is as expected, while 
the subleading divergence implies the 0(A 3 ) term of the /3-function: 

d ^ = -jk^ t) + TW m+0{>?) ■ (6 ' 10) 

Now we compute the momentum scale expansion non-perturbatively on (5.7). 
The 0(p°) case is evidently solved by 70 (pi, • • • , P4! t) = A(£), i.e. momentum 
independent, and gives the one-loop /^-function (5.8) as the flow equation. Following 
the prescription (4.5)-(4.8) we see that the 0(p 1 ) case generates, from the expansion 
of the propagator and I.R. cutoff, the non-analytic terms |pi + P2I, |pi + P3I and 
I Pi + P4I i-e- VS, VT and VU in Mandelstam variables. Therefore we try 

7(p 1 ,---,p 4 ;t)=A(t)+7 1 (t)(v / ^+v / T + v / t7) . (5.11) 

This gives, for example, (working always to order p 1 ): 

7(q> -q - Pi - P2, Pi, P2; t) = A + 2 7 i + ^71(1 + x) , 

where x = (pi + P2)-q/ V^- Substituting all expansions into (5.7) and evaluating the 
average we find that (5.11) is indeed the solution providing A(£) and 71 (t) satisfy: 
dX 3 



dt ' (An) 2 



(A + 2 7l ) 2 (5.12a) 



^ = -7i + ^3(A + 2 7 i)(A-[37r + 2]7i) , (5.126) 
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with boundary conditions A(0) = A and 71 (0) = 0. 

Integrating these equations numerically one finds that they focus in on the gaussian 
fixed point (A = 71 = 0) as required. This is illustrated in fig. 6. 
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Fig. 6. 71 plotted as a function of X(t) for bare couplings A = 20, 40, 100, 1000. 



For small enough A the solution for 71 is well approximated by the renormalization 
group improved perturbation expansion for this truncation, which is readily derived 
by using (5.12a) to solve (5.126) for 71 as a power series in A(t). One finds 



7i = 



1 a 2 + J-^A 3 + 



24tt 3 



96tt 5 ' 



+ 



1 



1536 7T 7 2304 7T 8 3456 tt 9 
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A 4 + 



(5.13) 



which substituted into (5.12a) gives the /3-function for this truncation: 

dX 3A 2 8 A 3 /32 16 \ A 4 

+ 5^ 7335-- • ( 5 - 14 ) 



dt (An) 2 n (4tt) 4 \n Sir 2 J (4tt) 

Several things should be noted about this solution. Firstly it is universal (i.e. 
independent of the cutoff). It may be shown that the non-perturbative corrections 
from (5.12) are all cutoff dependent. Secondly (5.13) and (5.14) are infinite series 
in the renormalised coupling constant. Therefore they incorporate non-trivial 
contributions to all loops. Thirdly, by comparison with (5.10) one sees that already 
at first order in the momentum scale expansion 99.1% of the two-loop contribution is 
included. In fact it may be shown that these series are asymptotic of the expected 
form, so they already capture much of the qualitative behaviour of perturbation 
theory to all orders. 

Now we briefly investigate the non-perturbative behaviour to test the convergence 
of the momentum scale expansion in this regime. Plotted in fig. 7 is a comparison 
between the zeroth order (5.8) and first order (5.12) results for A(t), given the 
extreme bare coupling Ao = 1000. 

We see that the difference is never large and falls to zero (as it must) as t — > oo. 
Since we are never interested in scales where cutoff effects dominate we quantify 
the difference in, what Liischer and Weisz refer to as, the "scaling region" where 
the computations are universal to a good approximation. For us this corresponds 
to the region X(t) < X s where, say, A s ~ 5 for Ao < 1000, since cutoff effects are then 
< 5%. (They reach ^5% at Ao = 1000 for A = 5, and are computed by comparing 
7i(t) with (5.13)). In this region all quantities may be expressed in terms of A(£), 
whose evolution is given by (5.14): 

16tt 2 8 , w . 
i = < ~ + 3AW + ^ lnA(<)+ ° <A) ' 

The integration constant too, the only remaining non-perturbative quantity, is 
readily computed numerically using the above formula. We find to first order in the 
momentum scale expansion = —2.98 for A = 1000, while too = —3.48 for A = 20, 
for example. These should be compared to the (zeroth order) one-loop result 
too = — 167r 2 /(3Ao). The effect on low energy is given by AA/A (16/t 2 /3)A Atoo- 
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Fig. 7. A comparison between zeroth A = \(°\t) and first order A = X^\t) 
approximations to X(t) in the momentum scale expansion, for bare coupling A(0) = 
A = 1000. Plotted are log 10 A^, log 10 A^ and AW/A (0) against t. 

Expressed as a percentage this is AA/A = aX where a = —5.6%, —1.7%, for Ao = 
1000, 20, respectively. 

Finally we briefly consider the more realistic case where a(p; t) is not set to zero, 
and (5.6) is therefore coupled to (5.1). From (5.11) one might expect a momentum 
scale expansion of the form cr(p;t) = o~o(t) +pai(t) + 0(p 2 ). If so, by dimensions 
one would expect o\ to be linearly divergent i.e. o\ ~ Ao/A. This would be a 
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disaster, as already explained at the beginning of sect. 4. In fact o\ is identically 
zero, since substituting (5.11) in (5.1) gives momentum dependence of the form 
< |q+ p| + |q— p| >q=i, yielding an expansion in p 2 . 

6. Summary. 

The main points of the paper are briefly recapitulated. 

The aim was to find a method of continuum calculation in realistic quantum field 
theories consisting of a sequence of better and better approximations, calculable 
without inhuman effort, and applicable even if there are no obviously identifiable 
small parameters to control the approximation. 

The Wilson renormalization group framework was chosen primarily because, unlike 
other frameworks, truncations are guaranteed to be at least perturbatively renor- 
malizable. 

Stated completely, the most promising method may be summarised as follows: Use 
the differential flow equation in A for the Legendre effective action 11 Ta[(/7 c ] defined 
by a theory with I.R. momentum cutoff A, and derived (formally) by ignoring the 
overall U.V. cutoff. Take the sharp cutoff limit in which the I.R. momentum cutoff is 
just a Heaviside function. ]?a[</? c ] may be thought of as equivalent to a Wilsonian 
effective action with effective U.V. cutoff A. The boundary condition to the flow 
equation is supplied at A = Ao by a rA [y c ], chosen to be local, which may be 
thought of as equivalent to a bare action. Finally, approximate the flow equation 
by truncation at some n-point 1PI (one-particle irreducible) greens function and 
expansion of each m-point greens function in momentum scale up to some maximum 
power. 

This prescription was arrived at as follows. The Wilson effective action £a [</?], 
where A is an effective U.V. cutoff, was shown to have a reinterpretation in terms 
of the generator of connected greens functions with I.R. cutoff A. As such it has 
an expansion in terms of 1PI greens functions, and thus is closely related to the 
Legendre effective action Ta[<p] in a theory with I.R. cutoff A. Unlike Sa[<p], Ta^] 

Here and later the effective actions are defined minus the classical kinetic term 
which incorporates the cutoff. 
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is insensitive to the precise form of the cutoff when the width e of the cutoff is 
small. It does appear to depend on an overall momentum cutoff Ao however, but 
we show that a tree- level reparametrization invariance allows us to change Ao. We 
use this invariance to interpret the flow equation for T\ [ip c ] in the simplifying limit 
A — > oo. This interpretation consists of defining the theory by a local 'bare' Ta [</? c ], 
the initial condition for the flow equation (from which a bare action and finite overall 
momentum cutoff can be reconstructed if desired). 

The 1PI greens functions, in the sharp cutoff limit e = 0, are non-analytic in the 
momenta around p = corresponding to non-local behaviour, however this appears 
merely to be a technical problem and not one of principle. On the other hand the 
smooth cutoff equations do not appear to have any useful approximation, the local 
derivative expansion having severe problems with convergence and sensitivity to the 
form of the cutoff. Straightforward numerical integration, expansion in ln(A /A), 
and expansion in A, are considered but are probably of limited use. 

Two simple model examples were considered to test the most promising method: the 
ladder (or rather cactus) approximation, and a model incorporating the first irrele- 
vant correction to the renormalized coupling, both for four dimensional A</? 4 theory. 
The cactus approximation is the simplest possible non-trivial non-perturbative 
approximation. It gives sensible qualitative results and even yields a maximum 
renormalized coupling in excellent quantitative agreement with previous precision 
work. For a certain range of bare couplings in the unstable symmetric phase it 
gives peculiar results, however it is easily seen that the approximation method does 
not converge there. The problem is avoided by shifting to the symmetric phase, 
or it may be resolved in other ways. In the truncation defined by setting the 
6-point vertex to zero, and further simplified by setting the self-energy to zero, 
the two-loop contribution was computed at zero external momentum using the 
momentum scale expansion, and seen to converge very rapidly, despite the fact 
that it is a numerical series with no obviously identifiable expansion parameter. 
Non-perturbatively the zeroth order in the momentum scale expansion coincides 
with the one-loop /3-function. First order in the momentum expansion gives a pair 
of flow equations whose (renormalization group improved) perturbation expansion 
yields an asymptotic power series for the /3-function, incorporating 99.1% of the 
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expected two-loop contribution. The non-perturbative corrections compared to the 
zeroth order result are small, giving confidence in a fuller use of these approximation 
methods for obtaining accurate results. 

A full demonstration that the proposed approximation method meets our require- 
ments of convergence and calculability awaits a proper calculation of some real 
non-perturbative problem. Research on using this to compute triviality bounds for 
the Higgs mass is underway. Since the method is tied to momentum cutoff, further 
work is needed before it can be applied to gauge theories. 

Note Added. 

Since this paper was submitted, a number of other relevant works have come to our 
attention. We here give a more complete comparison with earlier work. Of course 
there has been a vast amount of work on the exact renormalisation group in various 
guises; see for example ref.[l] for an early review and references. Continuum flow 
equations with smooth cutoff were derived in ref.[l] and are surely equivalent to 
those of Polchinski[2] (as indeed stated in ref.[2]). The sharp cutoff equations were 
derived by Wegner and Houghton[17], although the ambiguities we discussed in this 
paper were avoided by formulating the equations only for discrete momenta. It is 
this, presumably, that discouraged a more wide-spread use of these equations. The 
formal continuum limit for essentially these equations was derived in the appendix 
of ref.[18]. Weinberg noted that the equations were unpleasant, and by re-expressing 
the vertices as a tree expansion in some new vertices, derived flow equations for the 
new vertices in the sharp cutoff limit. The transformation in terms of trees was 
used in ref.[19] to establish equivalence between these equations and Polchinski's 
equations directly, for smooth cutoff also. Weinberg's equations coincide with the 
flow equations of the 1PI greens functions given in this paper (and refs.[5][8]), 
only the interpretation (in terms of greens functions as given here) is missing. 
Flow equations for the Legendre effective action (as given here and in refs.[5][8]) 
were derived in ref.[20]. I find it surprising that these equations appear not to 
have been used for approximations other than rederiving the e expansion [17] and 
constructing the simplest approximation where all momentum dependence 

is discarded [17] [6] -[9]. This latter equation was probably first explicitly written 
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down in ref.[21], the leading order large N case is given in ref.[17], see ref.[22] for 
a review. In this respect mention should also be made of Wilson's approximate 
recursion formulafl] which, like the p° approximation, also drops all momentum 
dependence, but here in an equation for halving the cutoff. As far as I know it gives 
similar results to p° approximation[7] [1], but in contrast to the p° approximation it 
cannot be regarded as the first term in a sequence of successive approximations [1]. It 
also differs in that the authors of ref. [1] prefer to consider smooth rather than sharp 
cutoffs. The reasons given are the difficulties inherent in using the Wegner-Houghton 
equation [17] and an understandable "philosophical" prejudice against the induced 
non-localities. As stated earlier in the paper, the transformation to 1PI vertices 
resolves the former's difficulties while, despite searching, we have found no deep 
problems with the latter. Also it should again be noted that differences between a 
sharp and smooth cutoff are in any case of no relevance to the p° approximation. 
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